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On the base of the quantum kinetic equation for density matrix in plasma at the stimulated 
bremsstrahlung of electrons on ions, the nonlinear absorption rate for high-intensity shortwave 
radiation in plasma has been obtained within relativistic quantum theory. Both classical Maxwellian 
and degenerate quantum plasma are considered for x-ray lasers of high intensities. Essentially 
different dependences of nonlinear absorption rate on polarization of strong laser radiation is stated. 
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I. INTRODUCTION 

At present, due to the rapid advance in free electron 
laser (FEL) technology EH became real the implemen¬ 
tation of ultrashort laser pulses in x-ray domain exceed¬ 
ing the intensities I0^° W/cm^. For such intensities the 
electron-wave field energy exchange over a x-ray wave¬ 
length is larger than the photon energy, and the laser- 
matter interaction has essentially multiphoton character 
Q- Thus, in recent decade the wide research field is 
opened up, where common nonlinear effects are extended 
to high energy transitions for various systems - from 
atoms [H, through molecules E3“E3l: ^o plasma and 
solid samples |14-[l^ . Among the fundamental processes 
of laser-plasma interaction, the inverse bremsstrahlung 
absorption of an intense laser field in plasma is one of 
the contemporary problems that have applications rang¬ 
ing from plasma diagnostics to thermonuclear reactions 
and generation of intense x-ray radiation. Note that in 
the field of intense x-ray laser, an electron may gain con¬ 
siderable energies absorbing even of a few quanta, which 
makes it as an effective mechanism for laser-plasma heat¬ 
ing. The theoretical description of this phenomenon in 
such superstrong radiation fields requires one to go be¬ 
yond the scope of common Quantum Electrodynamics 
-Feynman diagrams corresponding to the perturbation 
theory. 

With the advent of lasers many pioneering papers 
have been devoted to the theoretical investigation of the 
electron-ion scattering processes in gas or plasma in the 
presence of a laser field using nonrelativistic [20l - l3ll| as 
well as relativistic considerations. The appear¬ 

ance of superpower ultrashort laser pulses of relativis¬ 
tic intensities has triggered new interest in stimulated 
bremsstrahlung (SB) in relativistic domain, where inves¬ 
tigations were carried out in the Born (STI - K^ . eikonal 
, and generalized eikonal approximations over the 
scattering potential. Beyond this approximation for the 
infrared and optical lasers, in the multiphoton interac¬ 
tion regime, one can apply classical theory and the main 
approximation in the classical theory is low frequency or 
impact approximation . Low frequency approxi¬ 


mation have been generalized for relativistic case in Refs. 
[sM [3^ , where the effect of an intense EM wave on the 
dynamics of SB and nonlinear absorption of intense laser 
radiation by a monochromatic electron beam and by rel¬ 
ativistic plasma due to the SB have been carried out. 
For the infrared and optical lasers the quantum effects 
in the SB process is smeared out due to smallness of the 
photon energy. Meanwhile for a intense x-ray radiation 
the quantum nonlinear over the field effects will be con¬ 
siderable. Besides, relativistic quantum effects may be 
essential for plasmas of high densities [4l|, . Regard¬ 

ing the solid densities, absorption coefficient may reach 
considerably large values. As was shown in Ref. (H - 
where interaction of superstrong lasers with thin plasma 
targets of solid densities was investigated via particle-in¬ 
cell simulations, the inverse-bremsstrahlung absorption is 
dominant for electron densities above 10^^ cm“^. Thus, 
it is of interest to consider x-ray multiphoton absorption 
via inverse bremsstrahlung in ultradense classical as well 
as quantum plasmas. 

In the present paper the inverse-bremsstrahlung ab¬ 
sorption of an intense x-ray laser field in the dense classi¬ 
cal and quantum plasmas is considered in the relativistic- 
quantum regime considering a wave-field exactly, while a 
scattering potential of a plasma ions as a perturbation. 
The x-ray radiation power absorbed in plasma is investi¬ 
gated for a circularly and linearly polarized waves (CPW 
and LPW, respectivelly) arising from the second quan¬ 
tized consideration. The Liouville-von Neumann equa¬ 
tion for the density matrix is solved analytically for a 
grand canonical ensemble. With the help of this solu¬ 
tion we investigate the nonlinear inverse-bremsstrahlung 
absorption rate for Maxwellian as well as for degenerate 
quantum plasmas. It is shown that depending on the 
intensity of an incident x-ray, one can achieve the quite 
large absorption coefficients. Hence the considered mech¬ 
anism may serve as an efficient tool for ultrafast plasma 
heating. 

The organization of the paper is as follows. In section 
II, the relativistic quantum dynamics of SB is presented 
with analytical results for density matrix and inverse- 
bremsstrahlung absorption rate. In section III, we con¬ 
sider the problem numerically. Finally, conclusions are 
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given in Section IV. 


II. BASIC MODEL AND THEORY 

Let us consider the relativistic quantum theory of 
plasma nonlinear interaction with the arbitrary strong 
EM wave field by microscopic theory of electrons-ions in¬ 
teraction on the base of a single particle density matrix. 
The EM wave with the four-wave vector k = {lo /c, k) is 
described by the four-vector potential 

A'^ = (0,A), (I) 

where A is defined as: 

A(r) = Ao{ei coswT -|- e 2 gsina;T}; (2) 


, vor ck 

T = t - ; vq = —; eii'o = e2r'o = 6162 = 0, 

c w 

with the amplitude Aq, unit polarization vectors 61^2 
and ellipticity parameter g. The ions are assumed to be 
at rest and being either randomly or nonrandomly dis¬ 
tributed in plasma, the static potential field of which (for 
nucleus/ion -as a scattering center- the recoil momentum 
is neglected) is described by the scalar potential 

(x) = (:/9(r) ,0), 


where 


Ni 

(r-R-O- 
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Here ipi is the potential of a single ion situated at the po¬ 
sition Rj, and Ni is the number of ions in the interaction 
region. 

To investigate the quantum dynamics of SB we need 
the quantum kinetic equations for a single particle den¬ 
sity matrix, which can be derived arising from the sec¬ 
ond quantized formalism. As is known, the Dirac equa¬ 
tion allows the exact solution in the field of a plane EM 
wave (Volkov solution). Although the Volkov states are 
not stationary, as there are no real transitions in the 
monochromatic EM wave field (due to the violation of 
energy and momentum conservation laws), the state of 
an electron with a charge e and mass m in an EM wave 
field can be characterized by the quasimomentum 11 and 
polarization cr, and the particle state in the field ([ 1 ]) is 
given by the wave function: 




e(7fc)(7A) 

2c{kp) 


Uqjp) 


exp 



8c'^{pk) 


(1 — g^) sin(2a;r) 


(3) 


where V is the quantization volume, is the bispinor 
amplitude of a free Dirac particle with polarization cr, and 
n = {Eu/c, n) is the average four-kinetic momentum or 
’’quasimomentum” of the particle in the periodic field, 
which is determined via a free particle four-momentum 
p = {£/c,p) and relativistic invariant parameter of the 
wave intensity ^0 = eAp Imt? by the following equation 


2 2 

_ , m e , 9, ^9 


From this equation follows that 


m* = m [ I -b 


1+52 


1/2 


(4) 

(5) 


where m* is the effective mass of the particle in the 
monochromatic wave. The states © are normalized by 
the condition 

J = ^^^(5(n - n')d...'. 

Cast in the second quantization formalism, the Hamilto¬ 
nian is 


n = 


^+Ho^dr+nsb, 


( 6 ) 


where dr is the fermionic field operator, Hq is the one- 
particle Dirac Hamiltonian in the plane EM wave o, 
and the interaction Hamiltonian is 

^sb = lj (7) 

with the current density operator 

j = e$+7o7$. ( 8 ) 


Making Fourier transformation 


A(^)(x) = -^ [ (q)e-*‘i’’dq, 

j J 


the expression 0 will have a form 


Usb 


1 

c (27r)^ 


(q) e"*'i'’$dqdr. 


where 


(9) 


V{ci) 



R,) 


( 10 ) 


We pass to the Furry representation and write the 
Heisenberg field operator of the electron in the form of 
an expansion in the quasistationary Volkov states m 


eAp 
h [c{pk) 


(eipsinc 


ge 2 P coswr) 


§(r,t)=^ y'd$nan,<Te^^"‘d'n<7(i’,t), (11) 


X exp 
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where (i$n = Vd^II/(27rS)^. In Eq. (fTTI) we have ex¬ 
cluded the antiparticle operators, since contribution of 
electron-positron intermediate states will be negligible for 
considered intensities and photon energies = huj 
me?. The creation and annihilation operators, ^{t) 
and dn,cr(i), associated with positive energy solutions 
satisfy the anticommutation rules at equal times 

s (n - n') Sa.cr> , 

( 12 ) 

{®n,o-(^)> ®n',o-'(^ )}*=*' = {ciTi,(7{t), ajj',(j' {t )}t=t' = 0. 

(13) 

Taking into account Eqs. (fTTI) . ([5]), (O and the sec¬ 
ond quantized Hamiltonian can be expressed in the form 


+ie2g{Gs-i{a,l3,^) - Gs+i{a, P,^))} , ( 20 ) 

Bs = Gsia,p,^), (21) 

B2s = -^(1 + 9'^)Gs + “^(1 “ 5^) 

X {Gs-2{cx,/ 3,ip) + Gs+2ioi, I3,(p)), (22) 

and 

fiqs=n'-n-sfik (23) 

is the recoil momentum. The definition of the arguments 
a, P, If are: 


n^no + -Rsb (t). (14) 

The first term in Eq. m is the Hamiltonian of Volkov 
dressed electron field 


bAq 

he 


( flE 

V pk 


eiP'V 2 ( f2P 

p'k ) ^ \ pk 


e2p' 

p'k 


21 1/2 



Ho = 




d^uEua^ ^an,, 


(15) 


while the second term 

'Hsb (^) = ^ J J d‘hn'Afn',<T';n,<T (t) a^,^^,dn,r7 

go' 

( 16 ) 

is the interaction Hamiltonian describing the SB with 
amplitudes 


lbln'.cr';n,cr (t) 


- 


(17) 


8 ftc2 


(1-ff^) 




(25) 


tan(^ 


Q ( _ ^2P' \ 

^ \ pk p'k J 

( £1P _ ei P' ^ 

I pk p'k J 


(26) 


Thus, in order to develop microscopic relativistic quan¬ 
tum theory of the multiphoton inverse-bremsstrahlung 
absorption of ultrastrong shortwave laser radiation in 
plasma we need to solve the Liouville-von Neumann equa¬ 
tion for the density matrix p: 


M 


is) 


2cy^ EuEu’ 


ip') 


^='-[p,no + nsbit)], (27) 

with the initial condition 


/ eBiskep eepkBjA 
y 2c{kp') 2c{kp) j 


p{-oo) = pG- (28) 

Here pc is the density matrix of the grand canonical en¬ 
semble: 


e^(fceo)H2s - 
2c^{kp'){kp) 


Uaip)- 


(18) 


In Eq. (fTBl) the vector functions = (0, Bis) and 
scalar functions Hg, i? 2 s are expressed via the generalized 
Bessel functions Gs{a, P,(p): 


OO 

Gs{a,P,p)= ^ J2fe-.(a)Jfc(/3)e*(*-2'=)^ (19) 

k— — oo 



d^u {p — Eu)'> 


n, 


T^n,, 


(29) 

In Eq. (I^ni) Te is the electrons temperature in energy 
units, /i is the chemical potential, and H is the grand po¬ 
tential. Note that the initial one-particle density matrix 
in momentum space is 


Paic72(ni,n2,-oo) = Tt (pGa^^^^^dni,o 


Bis = ^ {ei (Gs-i(a,/3,v?) -b Gs+i(a,/3,/?)) 


= n (EnJ (ni-n2) 


V 


(30) 



















where 


X(5 {Eti' -EtiE sfvjj) (l - 


4 


n (Eyi^ ) 


1 


exp 


Eyii —M 

~^e 


+ l' 


(31) 


We consider Volkov dressed SB Hamiltonian 'Hgb (t) 
as a perturbation. Accordingly, we expand the density 
matrix as 


P = PG + '^^^- 


Then taking into account the relations 


®n' cr'^n.CTi Pg 


and 




j PG(^n',tT'^Tl,a 


[pGj'Ho] = 0, 


for we obtain 


in 


dt'^2 / / d^u'Mu'.,7>-,Ti,a (t') 


X e^(d-t){En,-En) -E^)'^ PG^^,^^^dn,a- 

(32) 

Now with the help of this solution one can calculate the 
desired physical characteristics of the SB process. In par¬ 
ticular, for the energy absorption rate by the electrons 
due to the inverse stimulated bremsstrahlung one can 
write 


^=Tr 

dt 


^1) 


dHsb (t) 
dt 


(33) 


It is more convenient to represent the rate of the inverse- 
bremsstrahlung absorption via the mean number of ab¬ 
sorbed photons by per electron, per unit time: 


dNje 1 d£ 
dt huNe dt ’ 


(34) 


where is the number of electrons in the interaction 
region. Taking into account decomposition 

(^1 - Tr (^ 001020304 ) 


= ni (1 - 712) i523(5i4, 

with the help of Eqs. (1321) . (l33l) . (l34l) . and (fTSl) for large 
t we obtain 


dNje 

dt 


dNje (s) 

^ It 

s—1 


(35) 


where the partial s-photon absorption rates are given by 
the formula 


dNje (s) Stts 


EnEw 


l^(q.)r 


B 


(s) 

n';n 


2 


xn{En')il-n{En)), (36) 


where 


B 


is) 

n':n 


e(pBiJa; e^ujB2s 
{kp)c 2c^{kp) 


|Bi«|'-Re(H2.B:)J|, 

(37a) 

and 5{x) is the Dirac delta function that expresses the 
energy conservation law in SB process. The obtained 
expression for the absorption rate is general and appli¬ 
cable to arbitrary polarization, frequency and intensity 
of the wave-field. This formula is applicable for a grand 
canonical ensemble and is always positive. With the help 
of Eqs. (1361) and one can calculate the nonlinear 
inverse-bremsstrahlung absorption rate for Maxwellian, 
as well as for degenerate quantum plasmas. 



III. NUMERICAL RESULTS AND DISCUSSION 


Eor the obtained absorption rate (1551) one need to con¬ 
cretize the ionic potential V (qs). For s single ion of 
charge number Za we will assume screening Coulomb po¬ 
tential with radius of screening as a function of the 
plasma temperature and density of electrons rig: 


Xe 



Thus, taking into account the plasma quasi-neutrality 
{ZaN^ = Ne), we have 


\V{qs)\^ = N, 


IGir^Zac'^ 


(38) 


Integrating in Eq. (l36l) over Eu' we will obtain the fol¬ 
lowing expression for partial absorption rates: 


dNrye (s) 2Zae‘^s 


dt 


TT^h^C* 


dEudfldfl' - 


iniin'i 


B 


is) 

n';n 


‘ m* c^-\-shLO 


-b 


X — e ^ ^ n {En — shu) (1 — n (En)), (39) 

where 

|n'| = ^\u\^ + h^s^k^-2^^^. 


dt 
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FIG. 1: (Color online) Envelope of partial rate of inverse- 
bremsstrahlung absorption via the mean number of absorbed 
photons by per electron, per unite time (in femtosecond”^) 
for CPW in Maxwellian plasma is shown for various wave 
intensities at = 1 keV, Za = 13, rie = lO^^cm”®, and 
Te = 100 eV. 



FIG. 2: (Color online) Same as figure 1 but for LPW. 


In general, the analytical integration over solid angles 
n, n' and energy is impossible, and one should make nu¬ 
merical integration. The latter for initially nonrelativistic 
plasma and at the photon energies hio > Te, is convenient 
to made introducing a dimensionless parameter 


Xo 


eEp 

ujy/rn}^^ 


(40) 


which is the ratio of the amplitude of the momentum 
transferred by the wave field to the momentum at the 
one-photon absorption. In Eq. (j40p the dimensionless 
parameter Eq = ujAq\/\ /c is the amplitude of the 
electric field strength. Hence the average intensity of the 
wave expressed via the parameter xo^ can be estimated 


FIG. 3: (Color online) Total rate of inverse-bremsstrahlung 
absorption for CPW in Maxwellian plasma versus the dimen¬ 
sionless parameter xo for various photon energies at = 
lO^^cm”®, and Te = 100 eV. 



Zo 


FIG. 4: (Color online) Same as figure 3 but for LPW. 


as 


=XoX 1-74 X 10^2 W cm-2 

The intensity strongly depends on the photon energy 
huj. At Xo 1) the multiphoton effects become essen¬ 
tial. Particularly, for x-ray photons with energies e-y = 
huj = 0.1 — 1 keV, multiphoton interaction regime can be 
achieved at the intensities ^ 10^® — 10^^ W/cm^. In 
the opposite limit xo ^ 1; the multiphoton effects are 
suppressed. 

For all calculations as a reference sample we take ions 
with Za = 13 (fully ionized Aluminum) and consider 
plasma of solid densities. In Fig. 1 and Fig. 2 the 
envelope of the partial rate of inverse-bremsstrahlung 
absorption for CPW and LPW in Maxwellian plasma 


huj 
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FIG. 5: (Color online) Total rates of the inverse- 

bremsstrahlung absorption scaled to Xo versus the parameter 
Xo for setup of Fig 3. 


is shown for various wave intensities at = 1 keV, 
He = and Tg = 100 eV (rig oc and 

Jig = (dTre^ng/Tg) ). As is seen from this figures, the 
multiphoton effects become essential with the increase of 
the wave intensity. 

To show the dependence of the inverse-bremsstrahlung 
absorption rate on the laser radiation intensity, in Fig. 3 
and Fig. 4 the total rate (1^ with (IMl) for CPW and 
LPW in Maxwellian plasma versus the parameter xo for 
various photon energies are shown. As is seen from these 
figures, the rate strictly depends on the wave polariza¬ 
tion, and for the large values of xo h exhibits a tenuous 
dependence on the wave intensity. 

To compare with the linear theory [^ . in Fig. 5 we 
plot scaled absorption rate Xo'^dN^g/dt versus xo- In the 
scope of the linear theory the scaled absorption rate does 
not depend on the wave intensity, while for the large val¬ 
ues of Xo it is suppressed with the increase of the wave 
intensity. To show the dependence of the considered pro¬ 
cess on the plasma temperature, in Fig. 6 we plot total 
rate of the inverse-bremsstrahlung absorption of circu¬ 
larly polarized laser radiation in Maxwellian plasma, as 
a function of the plasma temperature for various wave 
intensities at e.y = 1 keV, and Ug = 10^^cm“^. The 
similar picture holds for LPW. Here for the large values 
of Xo we have a weak dependence on the temperature, 
which is a result of the laser modified scattering of elec¬ 
trons irrespective of its’ initial state. We have also 
made calculations for a degenerate quantum plasma. In 
Fig. 7 and Fig. 8 the envelope of the partial rate of 
inverse-bremsstrahlung absorption for CPW and LPW 
in degenerate plasma with Fermi energy ^ ~ sp = 11.7 
eV (Aluminum) is shown for various wave intensities at 

= \ keV, and Tg = O.Ief (x'e = {Qi^e^ngj 
The total rate of inverse-bremsstrahlung absorption via 
the mean number of absorbed photons by per electron. 



Te [keV] 

FIG. 6: (Color online) Total rate of inverse-bremsstrahlung 
absorption of CPW in Maxwellian plasma, as a function of 
the plasma temperature is shown for various wave intensities 
at = 1 keV, Za = 13, rie = 10^®cm“®. 



Photon Number 

FIG. 7: (Color online) Envelope of partial rate of inverse- 
bremsstrahlung absorption for CPW in degenerate plasma 
with Fermi energy ep = 11.7 eV is shown for various wave 
intensities at = 1 keV, Za = 13, and Tg = O.Ief. 


per unite time for CPW and LPW in degenerate plasma 
versus the parameter xo at = HT eV and Tg = 0.1 
is shown in Fig. 9. As is seen from these figures, the 
rate strictly depends on the wave polarization, and for 
the large values of xo it is saturated. Note that our con¬ 
sideration is valid when the pulse duration r of an EM 
wave is restricted by the condition 

r < Vgff, (41) 

where Vgtf is the time scale during which the thermal- 
ization of the electrons energy in plasma occurs. In the 
presence of a laser field the electron-ion binary collisions 
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FIG. 8: (Color online) Same as figure 7 but for LPW. 



FIG. 9: (Color online) Total rate of inverse-bremsstrahlung 
absorption for CPW and LPW in degenerate plasma versus 
the parameter xo at sf = 11.7 eV, Za = 13, and Te = G.Ief- 


where Lcb is the Coulomb logarithm, and (v) is the mean 
values of electrons velocity in the laser field. For mod¬ 
erate intensities one can write (v) ~ fii® 

considered parameters we have — 10^"^ — 10^®s“^. 
Thus, the pulse duration should be r < 1 fs. The lat¬ 
ter is satisfied for x-ray sources. As is seen from Figs. 
3, 4, and 9, for the pulse durations r < 1 fs one can 
achieve an one absorbed x-ray photon by per electron, 
which means that in plasma of solid densities one can 
reach the plasma heating of high temperatures by x-ray 
laser with the intensity parameter ^ ^ 0.1 already at the 
one photon absorption process. 

IV. CONCLUSION 

Concluding, we have presented the microscopic 
relativistic quantum theory of multiphoton inverse- 
bremsstrahlung absorption of an intense shortwave laser 
radiation in the classical and quantum plasma. The 
Liouville-von Neumann equation for the density matrix 
has been solved analytically considering a wave-field ex¬ 
actly, while a scattering potential of the plasma ions as a 
perturbation. With the help of this solution we derived 
a relatively compact expression for the nonlinear inverse- 
bremsstrahlung absorption rate when electrons are rep¬ 
resented by the grand canonical ensemble. Numerical 
investigation of the obtained results for Maxwellian, as 
well as degenerate quantum plasmas at x-ray frequencies 
and large values of laser fields has been performed. The 
obtained results demonstrate that for the shortwave ra¬ 
diation the SB rate being practically independent of the 
plasma temperature is saturated with the increase of the 
wave intensity. The obtained results showed that in the 
x-ray domain of frequencies one can achieve an efficient 
absorption of powerful radiation specifically in plasma 
of solid densities and reach the plasma heating of high 
temperatures. 


take place with the effective frequency 


Veff ^ 


2TTZae^ne 
m? (v)^ 


Tcb; 
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